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Recently it has been suggested that many-body localization (MBL) can occur in translation- 
invariant systems, and candidate ID models have been proposed. We find that such models, in 
contrast to MBL systems with quenched disorder, typically exhibit much more severe finite-size 
effects due to the presence of two or more vastly different energy scales. In a finite system, this 
can easily create an artificial splitting of the density of states (DOS) into bands separated by large 
energy gaps. We argue that in order for such models to faithfully represent the physics of the 
thermodynamic limit, the ratio of the relevant coupling parameters must be larger than a certain 
cutoff that depends on system size, and should be chosen in such a way that various bands in the 
DOS of a given model overlap with one another. By setting the parameters in this way to minimize 
the finite-size effects, we then perform exact diagonalization studies of several translation-invariant 
MBL candidate models. Based on the variety of diagnostics, including entanglement properties and 
the behaviour of local observables, we find the systems exhibit thermal (ergodic), rather than MBL- 
like behaviour. Our results suggest that MBL in translation-invariant systems with two or more 
very different energy scales is less robust than perturbative arguments suggest, possibly pointing to 
the importance of non-perturbative effects which induce delocalization in the thermodynamic limit. 


PACS numbers: 73.43.Cd, 05.30.Jp, 37.10.Jk, 71.10.Fd 


I. INTRODUCTION 

One of the remarkable consequences of quantum me¬ 
chanics is the localization of a single particle moving in a 
disordered medium in one and two spatial dimensions [1]. 
This phenomenon, known as Anderson localization, re¬ 
lies on the presence of a static “quenched” disorder land¬ 
scape through which a particle can hop, causing its eigen¬ 
states to form localized wave packets instead of extended 
(Bloch) states in a usual crystal [see Fig. 1(a)]. This has 
dramatic consequences for dynamics and transport [2]. 

Recently it was found that localization can persist in 
interacting systems of many particles [3-6]. One of the 
simple examples of such systems is shown in Fig. 1(b); it 
consists of a ID lattice half-filled with particles that can 
hop as well as interact between nearest neighbour sites. 
Moreover, similar to the Anderson case, the system is 
in the background of an external disorder potential that 
is considered static and only enters the description as a 
random chemical potential on each site. Recent work [3— 
6] has shown that for sufficiently strong disorder such 
a system enters a “many-body localized” (MBL) phase 
where, similar to the Anderson case, transport is inhib¬ 
ited even at non-zero temperature. While in the Ander¬ 
son case the localization directly implies spatial locality 


of the wave function of typical eigenstates, in the inter¬ 
acting case each eigenstate |T) is a linear combination of 
Fock states |ni, n?, ..., n^), where n, is the number op¬ 
erator on site i. There are exponentially many possible 
|{rij}) and in a thermal (delocalized) system, |T) has a 
non-zero component on almost all of them. When dis¬ 
order is strong enough, it turns out that a typical |T) 
has non-zero projection on a far fewer {n,}). This is a 
signature of many-body localization, sometimes referred 
to as “localization in Fock space” because of the system’s 
inability to explore the entire many-body Hilbert space. 

A special structure of eigenstates IT) in the MBL phase 
generally is a result of an emergence of an extensive num¬ 
ber of locally conserved quantities [7-11]. The appear¬ 
ance of these local integrals of motion is responsible for 
the ergodicity breaking in the MBL phase, and leads to 
a characteristic entanglement structure [7, 12, 13] of its 
individual many-body eigenstates. In ergodic (thermal- 
izing) systems, the entanglement of typical eigenstates 
at high energy densities above the ground state obeys 
the so-called “volume law”: the von Neumann entropy 
of a large finite subsystem in such an eigenstate scales as 
the total number of degrees of freedom in the subsystem. 
This is simply a manifestation of |T) being a random 
linear combination of nearly all Fock states |{rij}). In 
contrast, the extensive number of local integrals of mo- 
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FIG. 1: (Color online) (a) Anderson localization of a single 
particle hopping on a lattice in the presence of external disor¬ 
der potential (blue circles). All eigenstates of such a particle 
in ID are exponentially localized in space, (b) Many-body lo¬ 
calization of particles hopping on a ID lattice and interacting 
between nearest neighbour sites. In this case also, the parti¬ 
cles move in the background of static external disorder poten¬ 
tial (blue circles). For sufficiently strong disorder, the wave 
function of typical eigenstates is localized (or weakly entan¬ 
gled) in the many-body Fock space, (c) A proposal for MBL 
in a translation-invariant system, similar to the ID Hubbard 
model. In this case, localization is expected to arise dynam¬ 
ically from strong on-site interaction between two species of 
particles with very different masses (A <C J <§C U), although 
no external randomness is imposed on the system. 


tion (which span the entire many-body Hilbert space) 
constrains the eigenstates in the MBL phase to have sig¬ 
nificantly “less” entanglement: their von Neumann en¬ 
tropy only scales as the area of the subsystem (there¬ 
fore, in ID the entanglement entropy is a constant). Fur¬ 
ther, dynamical probes such as the spreading of entangle¬ 
ment entropy following a global quench [14-16], as well 
as the time evolution of local observables [17-19], show 
that MBL phases have universal properties that distin¬ 
guish them from both the ergodic phase and the non¬ 
interacting Anderson insulator. 

One of the outstanding questions is whether the MBL 
phase necessarily requires disorder, i.e. could MBL- 
like physics and ergodicity breaking arise in translation- 
invariant systems, solely as a consequence of interac¬ 
tions? Recent works suggest that this intriguing phe¬ 
nomenon may indeed be possible [20-31]. A generic class 
of proposed translation-invariant MBL models involves 


two species of fermions a and b hopping on a ID lattice 
H = -~J^a\a i+1 - \^2b\b i+1 + h.c. 

i i 

+ TU(k- 0 a l a fc+<xft? bl+a '. (1) 

k,l 

Here J represents the hopping amplitude of fast (a) par¬ 
ticles and X J is the hopping of heavy (6) particles. 
The particles also interact via U(k — l). In the simplest 
case, cj = a' = 0, k = l and the interaction reduces to an 
on-site density-density term, thus the model is formally 
identical to the ID Hubbard model [30]. This model is 
schematically shown in Fig. 1(c) and can be viewed as 
a generalization of MBL models when the disorder is no 
longer static. That is, the particles that generate the dis¬ 
order are included in the system and allowed to undergo 
their own quantum dynamics and interact with the orig¬ 
inal particles of the system. Apart from the ID Hubbard 
model, we will also consider an alternative model with 
a = 1, a' = 0, k = l which can be visualized as light 
particles hopping on a lattice subject to the kinematic 
constraint depending on whether the hop extends across 
a heavy particle (“barrier”) or not [24]. Finally, we will 
also study a related one-component model which is the 
limiting case of the model a = 1, a' = 0, where light par¬ 
ticles have been integrated out, leading to an effective 
single species model with long-range interaction [29]. 

The expectation is that heavy particles in such models 
dynamically generate an effective disorder potential that 
localizes light particles. The possibility of such systems 
to display true [24, 29] or “partial” [30] MBL has recently 
been discussed. It was also argued that, regardless of 
whether such systems are truly MBL or not, they may 
break ergodicity in subtle ways that can only be detected 
by sophisticated probes such as the “post-measurement” 
entanglement of their eigenstates [23, 31]. 

In this paper we critically examine the mentioned fam¬ 
ily of translation-invariant candidate models for MBL. 
In Section II we consider models with two species of 
fermions (the ID Hubbard model and the model with 
particles and barriers), while Section III focuses on a re¬ 
lated single-particle model. Because the general theoret¬ 
ical understanding of such models is currently lacking, 
our approach combines insights from analytic solutions 
of small systems with exact diagonalization of larger fi¬ 
nite systems. Our conclusions are presented in Section IV 
and can be briefly summarized as follows. 

The main observation is that, in contrast to MBL 
models with quenched disorder, the translation-invariant 
models possess two or more drastically different energy 
scales (e.g., the hopping amplitudes of light and heavy 
particles). We argue that this leads to pronounced finite- 
size effects, and implicitly places constraints on the val¬ 
ues of parameters of the models (e.g., the minimal mass 
ratio of light and heavy particles). For example, for a 
given Hubbard chain of size L and fixed J, I/, there ex¬ 
ists a cutoff A C (L) such that only values of A > A C (L) 
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yield a faithful representation of the system in the ther¬ 
modynamic limit. For values A <C A C (L), the system’s 
DOS artifically splits into discrete bands separated by 
large energy gaps. The suppressed overlap between bands 
in DOS means that the effective Hilbert space is drasti¬ 
cally reduced. In this regime, a finite system indeed dis¬ 
plays MBL-like behaviour, as pointed out previously [30] . 
We argue that this feature, however, is a finite-size ef¬ 
fect that may not persist in the thermodynamic limit, as 
the broadening of each band is extensive in system size. 
Thus, the parameters of the model (e.g., A) have to be 
chosen carefully for a given finite system to avoid the 
mentioned DOS splitting. For such choices of A, we find 
that the typical behaviour of the system - diagnosed by 
the entanglement measures and thermalization of local 
observables [32-34] - appears to be ergodic rather than 
localized. Our results show that to demonstrate MBL in 
translation-invariant systems one needs to carefully ex¬ 
amine the finite sample’s DOS for the given hopping and 
interaction energy scales. More generally, the results im¬ 
ply that translation-invariant MBL is significantly more 
fragile than the disorder-driven one, and may generically 
exhibit a crossover to the thermal behaviour in the ther¬ 
modynamic limit. 



FIG. 2: (Color online) The DOS for the model (2) and 

separation into non-overlapping bands for L = 10, U = 10, 
J = 1. 


II. MODELS WITH TWO SPECIES OF 
PARTICLES 

We first consider a ID Hubbard model with two species 
of fermions a, b interacting on the same site, which is 
equivalent to the spin ladder model studied in Refs. [30, 


35]: 

H = ^ -Ja\a i+ 1 - Xb\b i+ i + h.c. + Up a i p \, (2) 

i 

where pf = crjai — a = a, b. The model is depicted in 
Fig. 1(c) and conserves the total number of particles of 
each species. We consider periodic boundary condition 
[i + L = i) and half filling for both a and b particles. We 
are interested in the regime A <C J -C 1/ where the system 
was conjectured to exhibit MBL behaviour [30]. We first 
discuss how the DOS in a finite system differs from ther¬ 
modynamic limit, and from this derive the conditions on 
the parameters of the model that minize the finite-size 
effects [Section II A]. In Section IIB we relate the be¬ 
haviour of entanglement and particle fluctuations to the 
discussed DOS and estimates derived in Section II A. 


A. —Finite system’s DOS and constraints on the 
parameters of the model 

The eigenstates of the model (2) have a simple struc¬ 
ture when particles are immobile (J = A = 0). In this 
case, a and b particles are randomly placed on lattice 
sites, and the energy is determined solely by U. The 
DOS in this case consists of massively degenerate bands 
that are separated by energy U. Next, let us turn on 
0 < J <C U while keeping A = 0. Now, a particles are 
Anderson-localized in the effective disorder profile gener¬ 
ated by the random positions of b particles. As a result, a 
band with a given classical interaction energy splits into 
many “mini-bands”. A typical mini-band still has exact 
degeneracy of at least 2 L, owing to the translation and 
reflection symmetry, and the typical energy difference be¬ 
tween neighbouring mini-bands is of the order A E(L), 
which decays exponentially in the thermodynamic limit. 
For J/U < 0.1 and the system sizes available numeri¬ 
cally (L < 10), the bands with different classical energies 
remain well-separated. In principle, this is a finite-size ef¬ 
fect that places the system far from the thermodynamic 
limit (effectively, because the system is too small to ab¬ 
sorb an energy of the order U). 

To illustrate the above, in Fig. 2 we show an exam¬ 
ple of the DOS for L = 10 site system and J/U = 0.1. 
Panel (a) is for the case A = 0.0005 where classical bands 
can be clearly identified. The bands are quasidegener¬ 
ate because of non-zero value of A; for A = 0 the bands 
indeed shrink to delta function peaks. In panel (b) we 
zoom in on one of the classical bands centered around 
e = —5. As stated above, we can now identify its fine 
structure corresponding to numerous mini-bands. Each 
of the mini-bands is at least 2 L degenerate, as we can 
see in panel (c) where we count the degeneracy of energy 
levels and plot the distribution of degeneracies (the lev¬ 
els are taken to be degenerate if their energy difference is 
less than 1CT 4 ). We also see that our estimate of typical 
degeneracy being 2 L is likely a conservative one because 
the distribution is peaked around the larger value which 
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in general is some multiple of 2 L [e.g., 4 L in Fig. 2(c)]. 
For comparison, in Fig. 2(d) we show the DOS for A = 1.0 
where the classical energy bands start to overlap. 
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FIG. 3: (Color online) The distribution of level spacings inside 
the central band of L = 10 system. We set U = 10 and A = 0, 
and assume periodic boundary conditions. The locations of 
the peaks are roughly consistent with the analytic estimate 
in Eq. (5). 

Let us provide an estimate for the critical A at which 
the DOS bands start to overlap. As each band contains 
an exponentially large number of states, turning on finite 
A can lead to thermalization within each band. We note 
that, since the bands are separated by an energy of order 
U , the hopping A typically couples states within the same 
band with an amplitude A J/U (since for many states 
moving a heavy particle changes the classical energy by 
U). If A is chosen very small, such that XJ/U is smaller 
than the mini-band spacing A E(L), then the effect of 
turning on the hopping is perturbative, and it cannot 
significantly modify the eigenstates. For instance, in our 
example in Fig. 2(a) we find the entanglement entropy of 
the eigenstates still obeys the area law. This is consistent 
with the system being in an MBL phase, but it is also 
likely that the finite system is then not reproducing the 
behaviour of the thermodynamic limit. Thus, to obtain 
a faithful representation of the thermodynamic limit in a 
large finite system of size L, we require A to be sufficiently 
large to hybridize several mini-bands 

A J/U > AE(L). (3) 

To estimate A E(L), consider the largest band with 
n(L) = (// 2 )(LL/ 4 J) s t a t es - At A = 0, a conservative 

estimate for the broadening of this band is (J 2 /U)'fL. 
We can easily convince ourselves that the broadening for 
A = 0 is given by ~ J 2 /U by exactly solving a small 
system, e.g. L = 4. Because of translation symmetry, we 
can restrict to a sector with the given momentum (e.g., 
k = 0). In this sector, there are 10 states whose exact 
eigenenergies are 

0, ±2J, ±s/U 2 + 4J 2 , ^{±U± \/U 2 + 16J 2 ), (4) 



thus the extremal eigenvalues are indeed shifted by ~ 
L(J 2 /U). In a large system, assuming normal distribu¬ 
tion for the modified energy values once we turn on J 
(while still keeping A = 0), we obtain the rough estimate 
for the broadening 


A E(L) 


2 J 2 L 3 / 2 
Un{L) ’ 


( 5 ) 


where we took into account the degeneracy of 2 L of a 
typical mini-band. For L = 10 this gives A E{L) ~ 
2.5 x 10 ~ 2 J 2 /U, which is of correct order but, as ex¬ 
pected, underestimates the exact result. Thus, the neces¬ 
sary (but possibly not sufficient) condition for the system 
to be able to capture the thermodynamic limit is 


A > A c « UAE(L )/ J « 0.025J. (6) 

In the following, we fix U = 10, J = 1. Signatures of 
MBL were found [30] for such parameters, when A was 
taken very small (10~ 3 —10 -2 ). Our estimate above gives 
A c ~ 0.025; choosing A as small as 10~ 3 would necessitate 
very large system sizes to make sure that finite-size effects 
are eliminated. 



FIG. 4: (Color online) Growth of entanglement entropy (left) 
and particle number fluctuations (right) in one half of the 
system as a function of time. System is L = 8 sites at half 
filling of species a and b. 


In Fig. 3 we test the derived estimate in Eq. (5) with 
the result obtained from exact diagonalization. We plot 
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the distribution of the level spacings between the mini¬ 
bands within the central (largest) band of the L = 10 
system (we set C/ = 10 and A = 0). This central band 
contains n(10) = 2520 states whose energy differences 
are plotted in Fig. 3. Level spacings between symmetry- 
related states have been removed, and the remaining 
ones are plotted on a log scale for easier comparison 
with our analytic estimate of A E(L) in Eq. (5). For 
the given choices of J, Eq. (5) predicts lnAEj =0 .ooi ~ 
— 19.8, lnAEj = o.oi ~ —15.2, lnAEj = o.i ~ —10.6 and 
lnAEj =10 ~ —6.0. These values are consistent with, 
though somewhat lower than the exact peaks of the dis¬ 
tribution in Fig. 3. 


B. Entanglement and particle fluctuations as a 
diagnostics of MBL 


Having derived the rough estimates for the parameter 
range in the Hubbard model, we proceed to compute the 
well-known diagnostics of many-body localization. One 
particularly useful probe is performing a global quench 
and examining the spreading of correlations after the sys¬ 
tem is initialized in a product state at time t = 0 and 
unitarily evolved with the Hamiltonian (2) for t. > 0. We 
will discuss the dynamical behaviour of the system un¬ 
dergoing a quench for the parameter values obtained in 
Section II A. 

In Fig. 4 we compute the evolution of entanglement en¬ 
tropy S ent under the global quench, as well as the particle 
number fluctuations (N%) — (1V CT ) 2 , er = a, 6, in one half of 
the system. Data is averaged over random initial product 
states. For an MBL system, we expect a slow, logarith¬ 
mic in time growth of S ent , and a much slower growth 
of particle fluctuations [15]. This behaviour reflects the 
slow dephasing and the suppression of transport, respec¬ 
tively. For the smallest value of A (roughly corresponding 
to Ref. [30]) we indeed find a signal that is reminiscent 
of MBL physics in the quenched disorder case [14-16]. 
The particle fluctuations, especially at long times, are 
significantly larger, which suggests the system becomes 
diffusive at very long times. However, for A > 0.1, where 
according to our estimate in Eq. (5) the system should be 
“closer” to the thermodynamic limit, we find signatures 
of thermal behaviour. For example, the growth of entan¬ 
glement is significantly faster and essentially featureless, 
with no indication of long time scales that characterize 
the transient regime. The saturation value of the entropy 
is high and approaches the thermal value independent of 
A for large A = 0.5,1, while particle number fluctuations 
grow at a similar rate for both species. These results 
suggest that as long as the DOS consists of overlapping 
mini-bands, the generic behaviour of the system is diffu¬ 
sive, and signatures of MBL only appear at intermediate 
time scales. 

We have reached similar conclusions by studying the 


model with particles and barriers proposed in Ref. [24]: 
H = ^2-J(a\a i+ 1 + h.c.){ 1 - b\bi) - Xb\b i+1 + h.c. (7) 

i 

This is the special case a = 1, a' = 0 of Eq. (1). In this 
model there is a kinematic constraint on the hopping of 
species a depending on whether the hop extends across a 
b particle or not. We have studied this model and found it 
displays similar phenomenology to the model in Eq. (2). 
These two-component models [Eqs. (2) and (7)], however, 
have the disadvantage of a rapid increase of Hilbert space 
dimension with system size L. Therefore, in the following 
section we focus on a single-component model [29] that 
is expected to capture similar physics, while at the same 
time allowing the numerics to reach bigger system sizes. 




FIG. 5: (Color online) The DOS for an L = 18 site chain 
at half filling described by the model (8). Decoupled bands 
are visible for small A = 0.001,0.01 (left), they begin to mix 
around A c = 0.1 and fully merge by A = 0.5 (right). 


III. A SINGLE-SPECIES MODEL 

Starting from the model (7), the strategy is to inte¬ 
grate out light particles a and obtain an effective single¬ 
component model for b particles only. This is possible 
under the approximation that there is exactly one light 
particle between the heavy ones [29]. In this case, we 
arrive at the following Hamiltonian: 

jj r ~ 1 

H = -A^6j5 i+ i + h.c. + rii'Y^ -jn i+r JJ(1 n i+fe ),(8) 

i r>0 k—1 

where rii = b] bi . Below we set U = 1 and (3 = 2 [29] . 

At A = 0, the DOS of the model (8) consists of several 
bands separated by gaps of order one. Each band con¬ 
tains configurations with the same number of occupied 
pairs of nearest-neighbour sites. The interaction energy 
is not exactly the same, as it depends on how the oc¬ 
cupied pairs of sites are arranged; thus, bands further 
split into mini-bands. In contrast to the Hubbard model, 
in this case the mini-bands can have higher degeneracy. 
Fig. 5 shows the DOS for L = 18 site chain at half fill¬ 
ing and different values of A. For small A = 0.001, 0.01 
(Fig. 5,left), the DOS consists of non-overlapping bands. 
By A = 0.1 these band start to mix, and for A = 0.5 the 
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DOS becomes contiguous. We therefore expect that crit¬ 
ical A c in this case lies in the interval (0.01;0.1). In the 
following we focus on three cases: A = 0.001, 0.1 and 0.5, 
and test whether local observables satisfy the “eigenstate 
thermalization hypothesis” [32-34, 36] (ETH), and study 
the scaling of entanglement in the eigenstates and its 
growth following a quench from an initial product state. 



0.2 0.21 0.22 0.23 0.24 0.25 

E/L 



0.2 0.21 0.22 . 0.23 0.24 0.25 

E/L 


FIG. 6: (Color online) Testing the ETH for the model (8) 
at half filling. Upper panels: Expectation value of the local 
operator O = nL/ 2 -\riL /2 in many-body eigenstates in the 
energy window 0.2 < E/L < 0.25, for several values of A. 
Bottom panels: the distribution of (O) for the same values of 
A. 

In Fig. 6 we test the ETH for the local operator 
O = nL/2-iriL/2 acting on the middle two sites in the 
chain. We compute the expectation value ( O) for the 
many-body eigenstates in the region 0.2 < E/L < 0.25 
near the middle of the band, for A = 0.001 and A = 0.1. 
Different choices of local operators O are possible, but 
yield qualitatively the same results. In the insets of Fig. 6 
we also show the distribution of ( O). For A = 0.001 the 
system is strongly non-ergodic, with a very broad distri¬ 
bution of ( O) [Fig. 6, top]. At A = 0.1 the system behaves 
according to the ETH prediction: the distribution of lo¬ 
cal observables becomes increasingly sharper with system 
size, and most of the eigenstates behave typically. 

Local observables therefore suggest that the model in 
Eq. (8) at A ~ 0.1 appears thermal rather than local¬ 


ized. We reach similar conclusions by considering nonlo¬ 
cal quantities such as entanglement entropy ^ent, shown 
in Fig. 7. Entanglement entropy S ent of typical many- 
body eigenstates is known to scale extensively with sys¬ 
tem size L in the ergodic case, and in contrast obeys an 
area-law scaling in the MBL phase [13]. In Fig. 7 we 
plot the entropy density S ent /L for all the many-body 
eigenstates within a given energy window in the middle 
of the band at A = 0.1. Note that we shift the data 
for different L by a certain amount E* dependent on L, 
so that the zero energy corresponds to the center of the 
band where states have the largest S ent . The characteris¬ 
tic concave shape of entropy is suggestive of thermal be¬ 
haviour; moreover, as we increase the system size, S en t/L 
appears to increase, which is inconsistent with localized 
behaviour where S en t/L should vanish. 



- 0.04 - 0.02 0 0.02 0.04 


(E - E*)/L 

FIG. 7: (Color online) Entanglement entropy density of the 
eigenstates in the middle of the band for various system sizes 
and A = 0.1 (left). Data for different L is shifted by E* 
so that zero energy corresponds to the center of the band 
where states have the largest S'ent • Average Sent in the window 
|E — E*\/L < 0.02 scales extensively for A > 0.1 (right). 

In Fig. 7(inset) we extract the average entropy (Sent) 
within the central window | E — E*\/L < 0.02, and per¬ 
form finite-size scaling for various A. Extremely small 
values of A such as 0.001 are consistent with area law 
that we expect in the MBL phase. However, for A > 0.1 
the average entropy within the central band clearly scales 
extensively with system size, suggesting that the system 
delocalizes. 

Finally in Fig. 8 we consider the dynamics of entan¬ 
glement entropy (left panels) and the particle number 
fluctuations (right panels) in the left half of the system 
for A = 0.001, A = 0.1 and A = 0.5. Particle num¬ 
ber fluctuations are defined as (N%) — ( Na ) 2 , where “A” 
stands for the left half of the system. It is instructive to 
consider both entropy and particle fluctuations because 
in the MBL phase with quenched disorder, the particle 
number fluctuations are strongly suppressed and depend 
weakly on system size due to the absence of transport in 
the system. On the other hand, the entropy still displays 
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FIG. 8: (Color online) Growth of entanglement entropy and particle number fluctuations in the left half of the system following 
the quench from a product state in the model (8). Data is averaged over random initial product states. Left panels show the 
entanglement entropy, right panels correspond to particle number fluctuations. Top to bottom: A = 0.001, A = 0.1 and A = 0.5. 


an unbounded logarithmic growth due to the clephasing 
dynamics in the MBL phase. In Fig. 8 finite systems of 
various sizes are prepared in an initial product state and 
evolved unitarily in time with the Hamiltonian (8). The 
data is averaged over random initial product states. Top 
row shows the case A = 0.001; for such a small value of 
A, the systems appears localized, although upon closer 
inspection we see that the long-time values of entropy 
and particle number fluctuations still grow slowly with 
L, indicating much stronger finite size effects than in the 
models with quenched disorder [16]. 

A more interesting case is A = 0.1 shown in the mid¬ 
dle row. If we look at Sent f° r a small system such as 
L = 8, we find some signatures of slow dynamics and 
“In t” growth of entropy prior to saturation, as expected 


of the MBL behaviour [15, 16]. However, we also no¬ 
tice that this slow growth of entropy is gradually washed 
away as we increase the system size to L = 18. For 
this large system, a much faster growth of entropy sets 
in, until the system starts to approach the saturation 
value of entropy. This faster growth is linear in time, 
as shown in Fig. 9, and is consistent with thermaliza- 
tion as we take L to infinity [37-39]. Similar conclusion 
can be drawn from particle number fluctuations [Fig. 8, 
right] which increase at long times, while they should be 
strongly suppressed in MBL phase with quenched disor¬ 
der [15]. From this, we conclude that A = 0.1 belongs 
to the transition regime between localized and thermal 
phases, where small system sizes L bias the systems to¬ 
wards the localized regime, but for sufficiently large L 
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FIG. 9: (Color online) Growth of entanglement entropy for 
A = 0.1 [middle row of Fig. 8] plotted on a linear scale. 

the system ultimately thermalizes. For comparison, we 
also show the case A = 0.5 [Fig. 8, bottom row] which is 
expected to be in the thermal phase. 

Finally one may expect in models such as (8) that the 
low density of particles may favor the localization. To 
address the dependence of our conclusions on the filling 
factor, we have also analyzed the case of 1/3 filling where 
we found qualitatively the same phenomenology as in the 
half-filling case: the entropy displayed faster than loga¬ 
rithmic growth in time upon the increase of system size, 
and the particle number fluctuations also continued to 
grow at long times as L approached the thermodynamic 
limit. 


IV. CONCLUSIONS 

We studied several candidate ID models proposed 
to exhibit many-body localization in the absence of 
quenched disorder. We emphasized that two or more very 
different energy scales in these models lead to pronounced 
finite-size effects, complicating the extrapolation of nu¬ 


merical results on small systems to the thermodynamic 
limit. More specifically, in a certain parameter regime 
(A -C 1), the DOS of small finite systems consists of 
bands separated by large energy gaps - a feature that 
generally does not persist in the thermodynamic limit. 
This artificial separation reduces the ability of the sys¬ 
tem to act as a heat bath for its parts. In this regime, we 
found, in agreement with previous works [29, 30], that 
the models exhibit phenomenology which is consistent 
with MBL. However, due to severe finite-size effects, it 
remains unclear whether such this type of MBL indeed 
persists in the thermodynamic limit. 

Further, we considered the parameter range A ~ 0.1 
where finite-size effects are weaker, and numerical simu¬ 
lations are expected to provide more reliable information 
about the thermodynamic limit. In this regime, using 
a variety of probes, we found thermal rather than MBL 
behaviour. Our results show that MBL in translation- 
invariant systems, if it exists, is much less robust than 
the perturbative arguments [24, 29] suggest (e.g., for 
the models studied here, perturbative arguments predict 
MBL, rather than thermal phase, at A ~ 0.1). Further 
work is needed to establish the mechanism of delocaliza¬ 
tion in these models, which might have non-perturbative 
origin [26]. 
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